Primordial black hole production due to preheating 
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During the preheating process at the end of inflation the ampHfication of field fluctuations can 
lead to the amplification of curvature perturbations. If the curvature perturbations on small scales 
are sufficiently large, primordial black holes (PBHs) will be overproduced. In this paper we study 
PBH production in the two-field preheating model with quadratic inflaton potential. We show that 
for many values of the inflaton mass m, and coupling g, small scale perturbations will be amplified 
sufficiently, before backreaction can shut off preheating, so that PBHs will be overproduced during 
the subsequent radiation dominated era. 
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I. INTRODUCTION 

There has recently been great interest in the reheat- 
ing process which occurs at the end of inflation. It has 
been found that in many inflation models reheating pro- 
ceeds initially via a period of broad parametric resonance, 
known as preheating During preheating the infla- 

ton field decays extremely rapidly producing an exponen- 
tially growing number of particles, until the backreaction 
of the particle production shuts off the parametric reso- 
nance with the particles produced subsequently thermal- 
ising and the universe becoming radiation dominated. 

Attention has been focused on two-field models with ei- 
ther a quadratic |^ or a quartic infiaton potential |^,^ . It 
has been claimed in Ref. Q that the amplification of cur- 
vature perturbations could be large with fluctuations be- 
ing driven non-linear, even on large 'supcr-Hubble' scales. 
This would have serious consequences for cosmology, de- 
stroying the standard picture of structure formation , 
in particular by violating the microwave background con- 
straint on the amplitude of the power spectrum on large 
scales, and leading to the overproduction of primordial 
black holes (PBHs) on small scales. The growth on large 
scales, relevant for structure formation, has been con- 
firmed for the quartic case |@-^, whilst it has recently 
been shown for the quadratic inflaton potential [p^|-p^, 
that in fact the amplification of the fluctuations on large 
scales is negligible. The formation of PBHs has not been 
studied for either model. For the quadratic potential 
the power spectrum of curvature perturbations has a k^- 
spectrum however, growing strongly towards small scales, 
so that PBH formation may still be problematic despite 
the amplification of large scale fluctuations being negli- 
gible. 

In this letter we present the first study of the formation 
of PBHs due to the resonant amplification of field fluctu- 
ations during preheating. The scales which we are inter- 



ested in pass outside the Hubble radius towards the end of 
inflation. The fluctuations are amplifled during preheat- 
ing and PBHs will be formed, when these scales re-enter 
the Hubble radius after reheating during the radiation 
dominated era, if the fluctuations are sufficiently large. 
PBHs may in principle also be formed on the slightly 
shorter scales which re-enter the Hubble radius during 
preheating, when the universe is dominated by the os- 
cillating scalar fields, however 3-D simulations along the 
lines of the ID simulations carried out in Ref. would 
be necessary to study this. 

II. CALCULATION OF THE POWER SPECTRUM 

In common with Refs. [§,0 we will study the two-field 
model with quadratic scalar field potential 

mx)-^mV + i5W, (1) 

where (j) is the inflaton and x the preheating fleld, m 
the mass of the inflaton and g the coupling of the in- 
fiaton to the preheat field. The cj) field, and hence the 
effective mass of the x field, = gcf), oscillates around 
zero with large amplitude resulting in efficient preheat- 
ing. The power spectrum of the curvature perturbation 
in the constant density gauge, C, produced in this model 
has recently been calculated by Liddle, Lyth, Malik and 
Wands [|2| (LLMW hereafter). The scales which we are 
interested in have k ^ 5$, where $ is the initial ampli- 
tude of the 0-ficld oscillations, and arc outside the Hub- 
ble radius at the end of inflation, so that we can use the 
results of their calculation. In this section we outline the 
relevant results of LLMW, for further details see Ref. . 

Strong parametric resonance, resulting in large ampli- 
fication of the initial quantum fiuctuation in the x field, 
occurs if 
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The power spectrum of a quantity x is defined as 
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where k = |k| is the comoving wavenumber, are the 
coefficients of the Fourier expansion of x, and the aver- 
age is over ensembles. The effect of preheating on the 
amphtude of the x field can be modelled as p|,p3 
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where At is the time elapsed since the end of inflation. 
The power spectrum of Sx at the end of inflation is de- 
noted by 'P^xlcnd- Floquet index fik is taken as 
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where fcmax is the scale which undergoes maximum ampli- 
fication and fcend is the comoving wavenumber of the scale 
which exits the Hubble radius at the end of infiation. For 
strong couphng » 1) Hk « /U = ln3/27r ft! 0.17 for all 
modes outside the Hubble radius at the end of infiation 
(k< 1). 

LLMW find that the evolution of the curvature pertur- 
bations is due to the non-adiabatic pressure perturbation 
i^Pnad, which is dominated by second-order fluctuations in 
X, g^'p'^^X^ (see also [0). The first-order contribution, 
which is of order g^^^X'^Xi is negligible since the back- 
ground field X is vanishingly small The power 
spectrum of ^nad resulting from the amplification of the 
non-adiabatic fiuctuations in the x field, is then given by 



A 



/(k, mAt) , 



(7) 



where 
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and 



I{K,mAt) = ^ dn' d6'e2(''»'+''«-')"'^*K'2sin6l, 

(9) 

where 9 is the angle between k and k' and Kcut is an 
ultraviolet cut-off, Kcut ~ '*max- LLMW found that on 
the scales relevant for large scale structure formation the 
change in C, due to preheating is negligible, however since 



'^Cnad smaller scale fluctuations with wavenumbers 

k ^ kcnd may become large enough for PBHs to be over- 
produced, before backreaction shuts off preheating. 

The duration of the first stage of preheating, mAtBR, 
during which backreaction is unimportant and fiuctua- 
tions are amplified exponentially, can be estimated nu- 
merically by finding the smallest solution of H 



lAt ft — In ^ 



(10) 



After this stage backreaction becomes important; energy 
is drained from the infiaton field, and transferred to the 
X field, diminishing the amplitude of the inflaton os- 
cillations. The growth in the preheating field changes 
the effective mass of the x field itself, thereby render- 
ing preheating less efficient. The growth in the x field 
also changes the oscillation frequency of the infiaton, nar- 
rowing the resonance band and eventually shutting the 
resonance down. The duration of this second stage of 
preheating is typically far shorter than that of the first 
stage of preheating. The solution of Eq. (p^, mAteR, 
thus provides us with a conservative estimate of the du- 
ration of the entire preheating process. 



III. BLACK HOLE ABUNDANCE CONSTRAINTS 

Primordial black holes may be formed in the the early 
universe via the collapse of sufficiently large density per- 
turbations ||l^ . Due to their observational consequences 
there are tight constraints on their abundance; typically 
less than 10~^° of the energy in the universe can go into 
PBHs at the time that they form ^j^. For a PBH 
to form during radiation domination the density con- 
trast, S, when the fiuctuation enters the horizon must 
exceed a critical value Jc, where numerical simulations 
find Sc ~ 0.7 iQ. The fraction of the energy of the uni- 
verse which goes into PBHs is found by integrating the 
probability distribution of the fluctuations, p{S), over the 
large, PBH forming, fluctuations. 



/3 = ^« / Pi6)dS, 

Ptot 



(11) 



and is controlled by the mass variance, (j{k), at horizon 
crossing, of the fluctuations: 
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The density perturbation distribution is close to gaus- 
sian; strictly speaking it is a chi-squared distribution 
with a large number of degrees of freedom, resulting from 
the convolution of the square of the gaussian distributed 
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FIG. 1. The mass variance at horizon crossing, (T{k), 
versus dimensionless wavenumber k/k^nd, where fccnd is 
the comoving Hubble radius at the end of inflation, for 
mAt — 0, 25, 50 and 75, for inflaton mass m = 10~^ and 
coupling g = 10""^. The dot-dashed line shows the threshold 
fthrcsh = 0.08 above which PBHs will be overproduced. 
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FIG. 2. The minimum preheating duration for which 
PBHs are overproduced, mAtpBH, versus the coupling g for 
q = lO*" and 10^. The dot-dashed lines show the time at which 
the first stage of preheating ends for the respective values of 
q. This shows that the fluctuations are amplifled sufficiently, 
before backreaction sets in, for PBHs to be overproduced sub- 
sequently. 



field fiuctuations over a wide range of scales.^ We find, 
using Eqs. (|ll|) and (|l^, that for a gaussian proba- 
bility distribution PBHs are overproduced, /3 > 10^"^", 
if cr(fc) > CTthrcsh = 0.08. For comparison for a first 
order clii-squared distribution the threshold would be 
Cthresh = 0.03. Wc usc (Tthrcsii = 0.08 as a conserva- 
tive estimate of the threshold for PBH overproduction. 
To calculate cr(fc) the power spectrum must be smoothed 
using a window function, W(kR), B: 
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The precise form of the window function is not important 
and we take it to be a gaussian. 



W{kR) = cxp 
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*This is significantly different from the case of single field 
models where difficulties arise in calculating /3 llSl]. In 
those models PBH production requires large field fluctuations, 
which are non-hnear and therefore may not be gaussian dis- 
tributed. Furthermore there is no definite prescription for 
relating the field and metric fluctuations if the field fluctu- 
ations are large. This is not the case in our model as the 
leading order metric perturbations arise from terms which 
are second-order in the linear field fluctuations. 



where R — 1/fc = 1/aH, is the scale of interest, i.e. the 
Hubble radius at the time the PBHs form. 

For mAt = 0, i.e. immediately after inflation ends be- 
fore preheating commences, /(k, 0) = k^^^ ^ 1, so that 
we can calculate (7{k) approximately analytically. Insert- 
ing Eq. (|^) into Eq. ( p^ ) and defining the dimensionless 
wavenumber k — k/k^aax we find 
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The integral in Eq. ([13) above can be evaluated as 
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so that the mass variance squared at the beginning of 
preheating is given by 
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\mAt=0 - ■■--1 1, 

which provides a check for our numerical results below. 
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IV. RESULTS 

We evaluate Eq. (Q) numerically using the values of the 
Hubble parameter H and the inflaton </> at the end of slow 
roll inflation: i/cnd = ?7i/3 and $ = 4>cnd — mpi/-\/127r. 

In Fig. |I| we plot the mass variance at horizon cross- 
ing, cr(fc), for several values of mAt, versus the dimen- 
sionless wavenumber k/kcnd, for the commonly used pa- 
rameter values m — 10~^mpi, g = 10~^; for these val- 
ues q « 6600. For mAt = we reproduce the ana- 
lytic expression found above, Eq. (|l^). We find that for 
mAt ~ 35 the mass variance a{k) exceeds (Tthrcsh — 0.08 
on small scales, leading to the overproduction of PBHs. 
Since a{k) is largest on small scales, PBHs are formed 
most easily immediately after preheating. In this pre- 
heating model k = aH oc a^^l"^ and the scale factor 
grows by at most 16 Q so that k changes by, at most, 
a factor of 0.25 during preheating. For definiteness we 
therefore evaluate a(k) on the scale fc* — 0.25fccnd- This 
approach is conservative; we over-estimate the Hubble ra- 
dius immediately after preheating and consequently also 
over-estimate the preheating duration neccesary to over- 
produce PBHs. 

In Fig. ^ we plot mAipBH, the minimum preheating 
duration for which the mass variance on small scales ex- 
ceeds the threshold for PBH overproduction, cr(fc*) > 
fthresh = 0.08, vcrsus the coupling for two values 
of the parameter q, q ^ 10*' and 10^. For comparison 
we also plot the estimate of the duration of preheating, 
mAteR, found from Eq. (|lO|). Similarly in Fig. ^we plot 
mAtpBH and mAieR as a function of g for m = 10~^mpi 
and 10~^mpi. The minimum preheating duration for 
which PBHs are overproduced, mAtpBH, depends most 
strongly on the strength of the coupling between the in- 
flaton and preheating fields, g, decreasing as the coupling 
increases. We find that, for all sets of parameters inves- 
tigated, the field fluctuations are amplified sufficiently 
during preheating that PBHs are overproduced during 
the subsequent radiation dominated era. The minimum 
preheating duration for which PBHs are overproduced 
is at least mAt 25 smaller than than the duration of 
the first stage of preheating, during which backreaction 
is unimportant and the fluctuations undergo exponential 
amplification. 
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FIG. 3. The minimum preheating duration for which 
PBHs are overproduced, mAtpBH, versus the coupling g for 
two values of the inflaton mass, m = lO'^mpi and 10~^mpi. 
The dot-dashed lines show the time at which the first stage 
of preheating ends for the respective masses. This shows that 
the fluctuations are amplified sufficiently, before backreaction 
sets in, for PBHs to be overproduced subsequently. 



fluctuations are amplified sufficiently during preheating, 
before backreaction is expected to shut-off the resonance, 
for PBHs to be overproduced during the subsequent ra- 
diation dominated era. We emphasize that, in the strong 
preheating regime, the overproduction occurs irrespec- 
tive of the strength of the coupling of the inflaton to the 
preheat field. This constitutes a serious problem for this 
preheating model. 
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V. DISCUSSION 



In this letter we have investigated primordial black 
hole formation in the two-field preheating model with 
a quadratic inflaton potential. Whilst the amplification 
of the field fluctuations on the scales relevant for struc- 
ture formation is negligible |10-I^, the power spectrum 



of curvature perturbations has a fc^-spectrum, i.e. the 
fluctuations are largest on small scales. We find that the 
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